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An ultra-short (about 30 fs) petawatt laser pulse focused with a wide focal spot (about 100 /mi) in 
a rarefied plasma (no ~ 10 17 cm -3 ) excites a nonlinear plasma wakefield which can accelerate injected 
electrons up to the GeV energy without any pulse channelling. In these conditions, propagation 
of the laser pulse with an over-critical power for relativistic self-focusing is almost the same as in 
vacuum. The nonlinear quasi-plane wake plasma wave, whose amplitude and phase velocity vary 
along the laser path, effectively traps and accelerates injected electrons with a wide range of initial 
energies. Electrons accelerated along two Rayleigh lengths (about eight centimeters) can gain an 
energy up to 1 GeV. In particular, the electrons trapped from quite a long (t& ~ 330 fs) non-resonant 
electron beamlet of 1 MeV particles eventually form a low emittance bunch with energies in the 
range 900 ± 50 MeV. All these conclusions follow from two-dimensional simulations performed in 
cylindrical geometry by means of the fully relativistic time-averaged particle code WAKE. 



PACS numbers: 52.35 Mw, 52.38 Hb, 52.38 Kd 
I. INTRODUCTION 



In the original, or " standard" , scheme of laser wake- 
field acceleration (LWFA) 0, , a laser pulse of duration 
smaller than a period of plasma oscillation r p = 2tt/lo p 
excites a wake electron plasma wave (wakefield) with a 
phase velocity close to the speed of light [here and else- 
where, uj p = (Airnoe 2 /to) 1 / 2 is the electron plasma fre- 
quency, tiq is the electron plasma density, to and — |e| are 
the electron mass at rest and charge] . When the electron 
density perturbation exceeds the background density, the 
accelerating electric field of the nonlinear wake can reach 
tens of GV/m, which is by three orders of magnitude 
higher than that can be achieved in the conventional ac- 
celerators without material breakdown. To excite such 
a wakefield, a short laser pulse of very high intensity is 
needed, that is, IqXq > 1.4 x 10 18 Wcm~ 2 /im 2 (where An 
is a laser wavelength). Reaching this range of intensities 
for lower laser energy in previous experiments Q made 
necessary a tight beam focusing. The laser focal spot 
with a radius ro ~ 10 /im gave a longitudinal extent of 
accelerating plasma (estimated as two Rayleigh lengths 
2Zr = 27rr 2 /A ) of the order of a few millimeters. Hence, 
the observed energy gain of electrons injected in the wake 
was moderate in spite of a high accelerating gradient. 
Also, very narrow plasma wake driven by a tightly fo- 
cused laser is ineffective for trapping and acceleration of 
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injected electrons: the scale of radial variation was con- 
siderably smaller than the axial wavelength, and most of 
the electrons were expelled by strong radial electric fields. 
Only few electrons were trapped and accelerated in the 
three-dimensional potential wells 3] . 

The petawatt (10 15 W) CPA lasers of new genera- 
tion, which are under construction now are capa- 
ble of delivering ultra-short pulses (tens of femtosec- 
onds) which can be focused with a relatively large focal 
spot (ro ~ 100 /im) yet having a relativistic intensity 
(Io > 10 18 W/cm 2 ) on axis. In this case, the effective 
acceleration distance (two Rayleigh lengths) is naturally 
elongated up to about ten centimeters, and the GeV en- 
ergy range can be achieved by using the standard LWFA 
technique without external optical guiding Q. In the 
scheme proposed, the plasma wavelength is smaller than 
the wake radial size. The large transverse extent of the 
wakefield structure reduces the detrimental effect of ra- 
dial forces on accelerated electrons. 

Anticipating this attractive chance to reach the GeV 
range of electron energy by using the standard LWFA, 
we analyze in this paper theoretically and numerically 
the details of laser pulse and wakefield evolution, and 
simulate the acceleration of an electron bunch in a wide 
range of parameters typical of those prospective experi- 
ments. The simulations have been carried out using the 
code WAKE 0, [|J . The laser parameters in the simula- 
tions are fixed: the laser wavelength Ao = 0.8 /im, the 
pulse energy 30 J, the pulse duration T£ = 30 fs, the 
radius of focal spot ro = 100 /im (hence, the maximum 
intensity in vacuum is Iq ss 6.4 x 10 18 W/cm 2 ). In the 
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simulations of electron acceleration the electron density 
of unperturbed plasma is n = 1.12 x 10 17 cm 3 , which 
gives the period of plasma oscillation r p w 330 fs and 
the plasma wavelength A p = 2ir/k p rs 100 /im (where 
kp = uj p /c). Then, the normalized pulse duration and 
width are u> p tl = 0.56 and k p ro = 6.28, respectively. 
The relativistic factor corresponding to the laser group 
velocity 7 g = ujo/lu p — 125. The laser power is exceeds 
by a factor of four the critical power for the relativistic 
self-focusing P cr = 16.2(w /w p ) 2 GW 0. 

The paper is organized as follows. Numerical study of 
nonlinear effects in propagation of ultra-short petawatt 
laser pulse through rarefied plasma is given in Sec. [H] 
The simulations supported by analytical considerations 
show that, for the parameters chosen, mutual cancella- 
tion of the relativistic and ponderomotive nonlinearities 
occurs in the pulse body. Thereby, only for very high in- 
tensities and/or electron densities (such that the critical 
power P cr is exceeded by more than a factor of four) the 
nonlinearity produces some effect on the pulse shape. 

Scction lTril addresses the effect of laser pulse transverse 
evolution on the phase velocity of the two-dimensional 
(2D) nonlinear wakefield. In our framework, the effect 
originates from the relativistic nonlinearity of the plasma 
wake. The amplitude of the wakefield is proportional to 
the laser intensity which grows (drops) while the pulse 
focuses (diverges). On the other hand, the relativistic 
reduction in frequency is proportional, in the weakly rel- 
ativistic approximation, to the wake amplitude squared. 
Therefore, as the pulse approaches the focal plane, the 
wake period stretches and its phase velocity drops be- 
low the pulse group velocity v g = c ho/wo. Beyond the 
focal plane the wake period shrinks and tends at infin- 
ity to the linear limit r p . At this stage, the wake phase 
accelerates and can become superluminous (i.e., exceeds 
not only the pulse group velocity in plasma but the vac- 
uum speed of light). A similar effect was discussed previ- 
ously for the case of longitudinally inhomogeneous plas- 
mas p|- Growth of the wake amplitude combined with 
the reduction in phase velocity provides conditions for 
highly efficient trapping of non-resonant injected elec- 
trons, 7 e o <C J g , in the stage of laser focusing [here and 
below, 7 e o = (1 — Vgo/c 2 ) -1 / 2 is the initial Lorentz factor 
of the injected electrons corresponding to the velocity i> e o 
at the point of injection]. 

Wakefield acceleration of resonant electron bunches is 
considered in Sec. IIVI Numerical simulations prove that 
the resonant conditions 7 e o ~ 7 9 or j c q w 7 p o (where 7 P o 
is a relativistic factor given by the local phase velocity of 
wake at the point of electron injection) are beneficial for 
the ultra-short electron bunches only that are loaded di- 
rectly in the accelerating and focusing quarter of the wake 
period (r p /4 w 80 fs for the parameters given). Taking 
a longer bunch leads to a considerable energy spread be- 
cause the resonant electrons injected at the bottom of the 
potential well stay there all the time and gain just a lit- 
tle energy. Injection of mono-energetic electron bunches 
with 7 e o = 7 9 = 125 or 7 e o = 42 in the second period of 



wakefield gives very similar output: after two Rayleigh 
lengths (about eight centimeters) the trapped electrons 
possess a broad energy spectrum (about 100% spread) 
with maxima at 0.5 GeV (for 7 e o = 125) or 0.75 GeV (for 
7 e o = 42) and a cutoff at approximately 1 GeV. Tuning 
the energy of injected electrons to the resonance with a 
given period of wakefield reduces the final rms emittance 
of the accelerated electron bunch but has no effect on the 
final energy spread. 

Acceleration of electron bunches injected with energies 
far below resonant is discussed in Sec. [V] Simulations 
show that the electrons with initial energies 5 MeV and 
1 MeV are not only effectively trapped and accelerated up 
to 1 GeV, but also reveal substantially less energy spread 
and lower root-mean-square (rms) emittance than in the 
resonant case discussed in Sec. IIVI The slow electrons 
(7e = 0.087 g and j eo = 0.0167 s ) loaded near the bot- 
tom of the 2D potential well slip into the accelerating 
and focusing phase of the wake period and get effectively 
accelerated. In particular, acceleration of the 1 MeV elec- 
trons along twice the Rayleigh length produces a group 
of electrons with the energy 900 ± 50 MeV and rms emit- 
tance less than 10 -4 mm mrad. 

Section IVII discusses and summarizes the results ob- 
tained. In Appendix lAl the analytical considerations are 
given of those nonlinear phenomena which have an effect 
on the propagation of ultra-short laser pulses in plas- 
mas. The features of the plasma wakefield generated by 
diffracting Gaussian short laser beam are considered in 
Appendix [5] 

II. SHORT LASER PULSE PROPAGATION IN 
TENUOUS PLASMA: COMPENSATION OF 
RELATIVISTIC AND PONDEROMOTIVE 
NONLINEARITIES 

Relativistic and ponderomotive nonlinearities of a 
short (ujpTL < 1) and broad (fc p ro > 1) over-critical 
laser pulse (P > P cr ) partly cancel each other The 
vacuum- like propagation of such pulse |10| . however, is 
not immune to the residual effect of the non-compensated 
relativistic nonlinearity (see Appendix lAl . We simulate 
the laser propagation, plasma wakefield excitation, and 
acceleration of injected test electrons by means of the 2D 
axially symmetric fully relativistic time-averaged particle 
code WAKE which is based on the quasi-static 

and extended paraxial |Tlj approximations. The normal- 
ized slowly varying amplitude (envelope) of laser vector 
potential a(£, r, z) is evaluated as a function of radial 
(r) and longitudinal (z) space variables and the retarded 
time £/c = t — z/c. In the simulation, the pulse prop- 
agates from the left to the right through a plasma slab 
centered at the vacuum focal plane z = 0. Simulation 
starts at a plane z = zq < and terminates at a symmet- 
ric plane z = \zq\. Initial laser beam with the parameters 
listed in Introduction is Gaussian [T^ . 

a = ao (r,e,Co)exp[-21n2(e-eo) 2 /(cTL) 2 + ^o] , (1) 



FIG. 1: The maximum intensity of laser pulse on the axis 
as a function of propagation distance. The vacuum focal 
plane is z — 0. The intensity is normalized to the maxi- 
mum intensity of laser pulse in the focal plane in vacuum, 
Io = 6.4 x 10 18 W/cm 2 . Solid line corresponds to focus- 
ing in vacuum; other lines correspond to focusing into plas- 
mas: dashed line — no/n c = 1.6 x 10~ 5 (P = P cr ), dash- 
dotted line — n /n c = 6.4 x 10~ 5 (P = 4P cr ), dotted line — 
n /n c = 1.28 x 1(T 4 (P = 8P cr ). 



where Co = -Zq/Z r , * = (r/r ) 2 Co/(l + Co)- arctan Co; 
a (r,(o) = a (l + C 2 )" 1/2 exp[-(r/r ) 2 /(l + Co 2 )] I £ = Co 
corresponds to the pulse center; and tl is the pulse 
full width at half-maximum (FWHM) in intensity. At 
the starting point, the laser pulse JU has a converg- 
ing (concave) phase front and focuses into plasma. The 
simulation proceeds from z = — Zr to z = Zr, where 
Zji w 4 cm for the parameters chosen. Normalized am- 
plitude of laser in the vacuum focal spot is ciq — 1.72. 

A few simulations with different plasma densities 
demonstrate the effect of nonlinearities on the pulse evo- 
lution. The maximum laser intensity on the axis (r = 0) 
is plotted in Fig.^versus propagation distance. The den- 
sity no = ngp = 2.8 x 10 16 cm -3 gives P — P cr (dashed 
line). At no = SnsF (or P = 8P cr , dotted line), the non- 
compensated relativistic self-focusing increases the inten- 
sity by a factor of 2.5 in the vicinity of the focal plane. 
However, for no = 4ns f (or P = 4P cr ), the on-axis laser 
intensity varies similarly to the vacuum case (dash-dotted 
line in Fig. [y. Therefore, under appropriate choice of 
parameters, propagation of a relativistically strong ultra- 
short overcritical laser pulse in a rarefied plasma is almost 
unaffected by the relativistic self-focusing, which could be 
a challenge for a longer pulse @ • 

Figure|5]demonstrates the distortion of the laser radial 
and temporal profiles due to the decompensation of rel- 
ativistic self- focusing. For P as 6P cr , the contour plots 
of normalized intensity a 2 (£, r) show that the pulse lead- 
ing part spreads in the course of propagation, while the 
trailing part shrinks. The effect manifests in full in the 
vicinity of the vacuum focal plane. The difference be- 
tween the speed of light in vacuum and the pulse group 
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FIG. 2: Contour plots of normalized pulse intensity at the 
three positions: z = — Zr (a), z = (b), z = Zr (c). At 
z — —Zr the pulse is Gaussian Q with the amplitude ao = 
1.72, radius k p ro = 7.84, and duration lo v tl = 0.68, which 
give P « QPcr- The solid, dashed and dotted-dashed lines in 
the plot (d) show the on-axis profile of intensity for the plots 
(a), (b), and (c). 

velocity in plasma brings about the gradually growing 
shift of the pulse center from its initial position in the 
co-moving variables. 

The analytical consideration given in a weakly rela- 
tivistic approximation in Appendix ^ predicts an am- 
plitude threshold ao c [see Eq. (|Allll ]. below which the 
laser pulse radially spreads according to the linear the- 
ory of diffraction for Gaussian beams For a Gaus- 
sian temporal profile, the critical amplitude reads ao c = 
(ujpTL,)- 1 x /41n2/[l + (fcpro/4) 2 ]. The laser amplitude 
ao = 1.72 taken in the simulation of Fig. [2] exceeds the 
critical amplitude ao c ~ 1.11, which explains the obvi- 
ous manifestation of the nonlinearity in Fig. [21 In order 
to reduce the adverse effect of decompensated nonlinear- 
ities, the wakefield evolution and electron acceleration 
will be further simulated in more rarefied plasmas. The 
parameters of the case P = 4P cr will be taken, that is, 
the plasma density no = 1.12 x 10 17 cm -3 that gives 
aoc ~ 1.6. The relativistic factor of the laser pulse then 
amounts to 7 9 = 125, and the normalized pulse duration 

tO LOpTL w 0.56. 



III. EXCITATION OF NONLINEAR PLASMA 
WAKEFIELD BY SHORT DIFFRACTING 
NEARLY GAUSSIAN LASER PULSE 

We study numerically and analytically the plasma 
wakefield evolution in the regime with laser nonlinearities 
mostly compensated and the laser beam close to Gaus- 
sian; the parameters are taken as suggested at the end of 
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FIG. 3: Temporal evolution of the radial profile of laser pulse 
intensity (in W/cm 2 ) (a) and electron density normalized to 
the background density no (b) in the vacuum focal plane z = 
0. 



section m Although the laser pulse is much shorter than 
a plasma period, the intensity on axis is high enough to 
produce in the wake an electron density perturbation of 
the order of the background plasma density. The wake 
becomes strongly nonlinear in the vicinity of focal plane, 
z = 0, where its amplitude is maximal. The radial and 
temporal profiles of the laser pulse and electron density at 
this plane are shown in Fig. [3J The wakefield is far from 
harmonic: the regions of density depression are much 
wider than the density humps. Moreover, the wake phase 
front is not plane and its curvature builds up with time. 
The relativistic shift of plasma frequency in the wake is 
proportional, to the lowest order, to the laser intensity 
squared, which radially varies and thus brings about the 
said curvature [liij . The radial phase variation, however, 
does not produce a radial wavebreaking [T3 | within at 
least six wake periods. 

The wake phase velocity is not constant along the laser 
path. This effect has a purely relativistic origin. As the 
pulse focuses at z < 0, its intensity grows, and so does 
the relativistic shift of plasma frequency. The wake pe- 
riod stretches, and the phase velocity drops below the 
group velocity of laser. Beyond the focal plane, z > 0, 
the pulse radially spreads, and the plasma period grad- 
ually shrinks thus tending back to the linear limit r p . 
At this stage, the wake phase velocity may exceed both 
laser pulse group velocity in plasma and the speed of light 
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FIG. 4: On-axis evolution of laser intensity (in W/cm 2 ) (a) 
and electron density normalized to the background density no 
(b) as a function of retarded time and laser pulse propagation 
length in the simulation of Fig. |3] The solid lines in the 
plot (b) are the tangents to the electron density crests. They 
characterize the local phase velocity of the wake at z = zo — 
— Zr. Corresponding relativistic factors in units of j g — 125 
are shown at the top line of the plot (b) . 



in vacuum. The luminous point, where the wake phase 
reaches the vacuum speed of light, does can exist on the 
path of a laser pulse in a longitudinally inhomogeneous 
plasma 8J. In our framework, it is due to the nonlinear 
frequency variation of plasma wake driven by the diffract- 
ing laser. Figure 0] shows the on-axis laser intensity (a) 
and the electron density (b) versus retarded time and dis- 
tance z in plasma. The plot (b) demonstrates the phase 
"deceleration" ("acceleration") at z < (z > 0) with the 
wake period stretching (contracting). At z > the wake 
phase is superluminous. 

An electron falls in resonance with the accelerating 
wakefield if its velocity v e o coincides at the point of in- 
jection with the local phase velocity of wake. Taking 
tangents to the electron density crests in Fig. Efb) helps 
to evaluate the resonant Lorentz factor j e o 3> 1 of the 
electron injected on axis at zq — —4 cm. The tangent 
equation, z w 2^7g , gives 7 e o ~ 42 s» 7 3 /3 for the sec- 
ond, and 7 e o w 21 = 7g/5 for the fifth wake period. 
Weakly nonlinear theory of the wakefield excitation by a 
diffracting short Gaussian laser pulse (see Appendix |B|) 
expresses tThe wake Lorentz factor as a function of nor- 
malized pulse position £ = z/Zr and time delay £/c [see 
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FIG. 5: Acceleration of the resonant electron bunch (7 e o = 
7 9 = 125). The space distribution (a), momentum distri- 
bution (b), and energy versus the injection phase (c) are 
shown for test electrons crossing the extraction plane z = Zr ■ 
Each dot corresponds to a numerical electron. The density of 
dots characterizes the number of test particles per elemen- 
tary volume kpdrd^ in the plot (a) and per elementary vol- 
ume of momentum space (dp r dp z ) / (mc) 2 in the plot (b). The 
grayscale density in the plot (a) is proportional to the normal- 
ized wakefield potential (ip), the lightest gray corresponds to 
(V')max = 0.3, the darkest gray — to (i/Omm = —0.24. Dashed 
lines in the plot (a) show the boundaries of the bunch at the 
injection point. 

Eq. I|B5|) ]. which gives 7 e o by a factor of two higher than 
inferred from the graphical estimates of Fig. 0Jb) . 



IV. ACCELERATION OF RESONANT 
ELECTRON BUNCHES 

The test electron bunch is injected in the laser wake at 
zq = —Zr. At this point, a Monte-Carlo generator cre- 
ates mono-energetic particles with the energy mc 2 7o e and 
uniform distribution over time interval Tf,. The bunch has 
a transverse momentum spread, which gives a nonzero 
initial angular divergence (cv 2 ) 1 / 2 = a e , and a nonzero 
initial rms emittance e± = (l/2)[(r 2 _)((pj_/p) 2 ) — ((rj_ ■ 
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FIG. 6: Energy spectrum (a) and emittance (b) at the ex- 
traction plane z — Zr for electron bunches injected with 
7 e o = 125 (thick lines) and 7 e o = 42 (thin lines) into the 
second wake period. 

Px)/p) 2 ] 1//2 (Ref. \T3 ). Radial distribution of electron 
density in the bunch is Gaussian with an rms radius 
a ~ tq. The accelerated electrons are extracted from 
the wake at z = Zr. 

First, the conventional resonant condition 7 e o = "f g is 
considered. A bunch of 5000 test electrons with j e Q = 
125, zero angular divergence, and k p a = k p r /^/2 = 4.5 
is accelerated in the second wake period (6.5 < k p ^ < 13, 
which corresponds to the bunch duration r& ps 330 fs). 
Distributions of test electrons in coordinate and mo- 
mentum space at z = Zr , and final energy versus 
injection phase are shown in Fig. [3] The grayscale 
background in the plot (a) is the normalized potential 
(-0) = [e/(mc 2 )](A z — <£>), where (• • • ) means the averag- 
ing over the laser period 2tt/ujo, (A z ) and ($) are vec- 
tor and scalar potentials associated with low-frequency 
wakefields. This potential determines the low-frequency 
electric and magnetic fields in plasma, and, hence, the 
forces acting on test electrons. The plot (a) shows that 
accelerated particles are collected in the focusing and ac- 
celerating quarter of the wakefield period. Quite a few 
electrons stay near the bottom of the potential well with 
the energy gain, as is testified by the plot (c), almost 
negligible. The transverse spread of the bunch is reduced 
roughly twice if compared with the initial value. 

Figure [S] shows the final energy and emittance of the 
test electrons collected by equally distributed channels 
of an electron spectrometer. We compare them for the 
two different resonant condition fulfilled at the injection 
plane. The thick lines correspond to 7 e o = 7 9 = 125 
(resonance with the laser pulse) and the thin ones to 
7 e o = 7po ~ 42. This value of "f e Q corresponds to the 
resonance with the second period of the plasma wake 
when the injection point is on the axis. In both cases 
the energy spectrum looks as a shoulder with not too 
pronounced maxima near 0.5 GeV for j e o — 125 and 
0.75 GeV for 7 e0 = 42. There is no explicit energy gap 
separating the accelerated electrons from those which are 
non-accelerated. In the case of j c q — 42 the final emit- 
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tance is smaller than that for 7 e o = 125 [see the plot 
(b)]. The energy cutoff appears to be independent of the 
injected electron energy and is close to 1 GeV. 

The energy limit of 1 GeV can be exceeded if the elec- 
trons are injected earlier and extracted from the wake 
later. Additional runs show that elongating the acceler- 
ation length by 50% (-1.5Z R < z < 1.5Z R ) gives a 15% 
increase in the maximum of final energy. As the highest 
acceleration gradients are achieved near the laser focus, 
increasing the acceleration length beyond 2Z R does not 
lead to a substantial growth of electron energy. 

Analysis of Figs. [5] and leads to the following conclu- 
sion: injection of a quite long (of the order of or longer 
than a wake period) resonant electron bunch results in 
a large energy spread of trapped and accelerated elec- 
trons no matter whether the injected bunch was resonant 
with the laser pulse or with a given period of the wake. 
It is obvious that the resonant conditions "f e o = j g or 
7 e o = 7 P o serve well only for the electrons injected in just 
one quarter (focusing and accelerating) of the wake pe- 
riod. Sample simulation made for a short electron bunch 
(11.5 < kpt; < 12.5) with 7 e o = 125 and zero initial emit- 
tance gave the mean output 0.85 GeV per electron with 
about 30% energy spread. 



V. ACCELERATION OF NON-RESONANT 
ELECTRON BUNCHES 

Now studying non-resonant injected electron bunches, 
we bear in mind that the initial energy of electrons is es- 
sentially smaller than that of the resonant ones. Figure[7| 
presents the same plots as Fig. 5 but for injected electrons 
with 7 e0 = 10 < 7 9 (p0)< 

Comparison of Figs. and reveals certain benefits 
of using the electrons non-resonant at injection. First, 
the accelerated bunch of initially non-resonant particles 
is more compact in the radial direction [compare plots (a) 
of the two figures] , and the radial momentum spread is 
lower [compare plots (b)]. Fig.[7{a) shows no electrons at 
the bottom of the potential well, « 10, hence, all the 
trapped non-resonant electrons are accelerated, and the 
energy gap finally appears between the trapped and non- 
trapped particles [see Fig.UJc)] . Opposite to the resonant 
case, the non-resonant electrons, loaded into either accel- 
erating or decelerating focusing phases, gain almost the 
same energy. At the point of injection, the wake phase 
outruns the non-resonant particles, so they slip out of the 
disadvantageous phase towards the focusing and acceler- 
ating one. Otherwise, the resonant electrons are initially 
at rest in the decelerating phase, and their slippage takes 
more time. So, their final energy proves to be lower. 

Figure [S] shows final energy spectrum and emittance 
of initially non-resonant test electrons with j e Q = 10, in- 
jected at zq = —Z R with zero initial omittance in the 
second (6.5 < fc p £ < 13, simulation of Fig. UJ and sixth 
(32.1 < kpt; < 38.4) wake periods. The initial rms radius 
of the bunch is the same as in Figs.[S]andU| The electrons 
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FIG. 7: Acceleration of the non-resonant electron bunch. All 
parameters being the same as in Fig. except of initial rela- 
tivisitc factor of electrons 7 e o = 10 <C y g . 



accelerated in the second wake period are separated from 
the bulk of low energetic particles by a 150 MeV wide gap 
[see also Fig.Ufc)]. Besides, efficiency of acceleration re- 
duces when the time delay between the laser pulse and 
electron bunch grows. The spectrum of electrons acceler- 
ated in the second period has the maximum at 0.8 GeV 
and a sharp cutoff at 0.95 GeV. Acceleration in the sixth 
period gives a plateau-like energy spectrum which rises 
steadily up to 0.65 GeV and drops at 0.85 GeV. There- 
fore, the injection time lag should be taken as small as 
possible in order to reduce the adverse effect of the wake 
phase velocity variation. As compared with the resonant 
case [Fig. EJb)], final emittance of the bunch is typically 
lower in the non-resonant case (by a factor of three for 
the spectrometer bins beyond 500 MeV). 

Technological limitations of monoenergetic electron in- 
jectors lead inevitably to using in the experiment long 
(t& ^> t p ) and wide (k p a ^> 27r) electron bunches with 
nonzero rms emittance. The results of sample modelling 
of this situation are presented in Fig.|5| The energy spec- 
trum (a) and rms emittance (b) of accelerated electrons 
are shown at the extraction point z = Z R . The non- 
resonant (7 e o = 10) electron beam injected at zq = —Z R 
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FIG. 8: Final energy spectrum (a) and emittance (b) of elec- 
trons injected in different wake periods. The electron are 
injected with 7^0 = 10. Thick and thin lines correspond to 
the injection to the second (data from simulation of Fig. 
and sixth period, respectively. 
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FIG. 9: Energy spectrum (a) and rms emittance (b) of ac- 
celerated electrons for non-resonant (7 e o = 10) wide (k v a = 
13.5) and long (6.5 < k p £ < 25.4) electron bunch with zero 
initial emittance (thick lines) and with an initial rms emit- 
tance £j_ = 2.145 mm mrad (thin lines). 
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FIG. 10: Acceleration of low-energy electron bunch (7 e o = 2). 
Plots (a)-(c) show the same quantities as in Figs. and 
The electron energy spectrum is shown in plot (d). Injection 
of highly non-resonant bunch provides almost monoenergetic 
acceleration of trapped electrons. 



eV], which lead to almost monoenergetic accelera- 

Eg. In 



1 MeV). 

tion [lg. In Fig. EH the energy spectrum of accelerated 
electrons with initial relativistic factor 7 e o = 2 and radial 
spread k p a — 2 is shown while the other parameters are 
taken the same as of Figs. |5] and The total number 
of injected particles is 2500, of which roughly one half is 
trapped and accelerated up to 900 MeV with 10% energy 
spread. The final rms emittance of accelerated bunch is 
lower than 10~ 4 mm mrad. 



covers three consecutive wake periods, 6.5 < k p £ < 25.4. 
The transverse size of the beam k p a = 13.5 is by a factor 
of three larger than that of Figs. and and the num- 
ber of particles in the bunch is 20000. The thick lines 
correspond to the bunch with zero rms emittance at in- 
jection, and the thin lines to the case with initial rms 
emittance e± = 2.145 mm mrad. Variation of the initial 
rms emittance has a negligible effect on the energy spec- 
trum which has a sharp maximum at 0.75 GeV and cut- 
off at 0.95 GeV. Only emittance of low-energy electrons 
(spectrometer bins at E < 500 MeV) is increased due to 
nonzero rms emittance of the injected beam. The emit- 
tance of the higher-energy electrons is unaffected. Thus, 
the simulations show that elongating the electron bunch 
still preserves all the benefits of using the non-resonant 
electrons and does not create any challenge for the ex- 
perimental implementation of the proposed scheme. 

Effective trapping and acceleration of electrons also oc- 
curs for bunches injected with very low energy {Ei n ~ 



VI. CONCLUSION 

Construction of the petawatt ultra-short-pulse lasers 
will create an opportunity for standard LWFA of elec- 
trons up to GeV energies. In quite tenuous plasmas 
(uip < t7 ) the relativistic and ponderomotive nonlin- 
earities of a not very tightly focused (k p r > 2n) laser 
pulse cancel each other |9j. In this regime, an overcrit- 
ical laser pulse propagates like in vacuum. Increasing 
the laser focal spot up to 100 /im in radius elongates the 
laser-plasma interaction length (estimated as twice the 
Rayleigh length) up to 8 cm without any external opti- 
cal guiding. Such not very tightly focused laser pulse, 
however, has enough energy to drive a nonlinear plasma 
wakefield along this distance, which provides a control- 
lable acceleration of externally injected electrons up to 1 
GeV. The nonlinear features of quasi-plane wake plasma 
wave facilitate trapping, focusing, and acceleration of the 
electrons from the injected bunch. The converging laser 
pulse drives the plasma wake whose period stretches (due 
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to relativistic decrease in plasma frequency) as the pulse 
approaches the focal plane. At this stage, the wake phase 
velocity drops below the pulse group velocity. Growth 
of the wake amplitude combined with the decrease in 
the phase velocity provides the efficient trapping of low- 
energy (non-resonant) electrons. Beyond the focal plane, 
the laser pulse diverges, and the wake phase may become 
superluminous, which makes for additional longitudinal 
compression of the electron bunch in the focusing and 
accelerating quarter of the wakefield period. These fea- 
tures of wakefield structure reduce the energy spread and 
omittance of the electron bunch injected out of resonance 
with the wake, 7 e o <C 7 9 (po): an d, m the limit of very low 
initial energy (1 MeV), provide almost monoenergetic ac- 
celeration. 
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APPENDIX A: NONLINEAR EFFECTS IN 
PROPAGATION OF ULTRA-SHORT LASER 
PULSES 

Self-consistent evolution of axi-symmetric laser pulse 
and perturbations of electron plasma density can be de- 
scribed in the weakly nonlinear quasi-static approxima- 
tion by the set of equations |l7j. 



2ik — + Aj_ 
dz 



= k 



1 



•V--|o| 
A j 



" V 4 
1 / d 2 



(Ala) 
(Alb) 



where a — eE/(mujQc) is the normalized amplitude 
of the laser electric field, N = (n e — no) /no is the 
normalized electron density perturbation, and Aj_ = 
r~ 1 d/dr(rd/dr) is the radial part of the Laplace opera- 
tor. 

In the limit of a short and wide pulse whose length 
ctl is smaller but the radius ro is larger than A;" 1 , the 
main terms in Eq. I jAlbj l are those containing the sec- 
ond derivatives with respect to £. In this approxima- 
tion, the nonlinear terms in the right-hand side (RHS) 
of Eq. (|Ala(l cancel each other, and the pulse propagates 
as in vacuum |9(. The small corrections to the index of 
refraction which are due to the finite pulse length and 
proportional to {lo p tl) 2 were considered in Ref. 0- We 



consider the effect of small longitudinal and transverse 
nonlinearities on the evolution of a short laser pulse and 
assume the transverse pulse shape be Gaussian in every 
cross-section, 



MO 

a=—— — exp 

/(£>*) 



,2 /2 



(A2) 



Here, /(£, z) is the dimensionless focal spot size (pulse 
width) which equals unity in the focal plane z = [i.e., 
/(£,0) = 1]. The function A (0 gives an initial ampli- 
tude profile on axis, the initial focal spot radius is r , 
the quantities tp(£, z) and /?(£, z) give the on-axis values 
of phase and curvature of the laser phase front, respec- 
tively. Substituting Eq. (|A"2l) into Eqs. {£H § ives tlie 
equation for the pulse width [l8| . 



dz 2 



7 2 



4 u) Q Z R 



(A3) 
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f*(£>,z) V ' k 2 p r 2 f 2 (^,z) 



where z is the coordinate of the pulse center moving from 
the left to the right; the transverse Laplace operator in 
Eq. jug gives the second term in the brackets in the in- 
tegrand. In the short-pulse case, Aq(£) is nonzero within 
an interval |£| k^ 1 . Integrating by parts the expres- 
sion in the RHS of Eq. (|A3(I and taking account of only 
the main term gives 



dz 2 



7 2 



4 lo Z r 



(A4) 



16 



The laser pulse of small amplitude propagates in plasma 
as in vacuum. In order to evaluate the threshold ampli- 
tude aoc above which the effect of nonlinearities might 
occur, we consider Eq. (|A4|) near the focal plane. At the 
focal plane z = the pulse width / is constant while 
alters within a relatively short interval |£| <C k^ 1 . 
We assume that f(£,z) varies with £ slower than j4o(£) 
within some segment of the laser path. Equation (|A4|) 
then reduces to 



(A5) 



where £ = z/Zr is the normalized propagation distance, 
a = (fc p r ) 2 /8 is the normalized initial squared spot size 



no = k 2 P f ^ [ dCAKc 

J — oo J — oo 



(A6) 



The condition /(e, 0) = 1 gives the integral of Eq. I|A5(1 . 

2 



d() \ f 2 



I -nl - / ( 1 + JL 



= 0. (A7) 
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Initially, the evolution of a small spot size perturbation 
5f = f—l,\Sf\<l, obey the equation following from 
Eq. f£7L 



(d6f/dq) 2 -2Sf(B-l) = 0, 



(A8) 



where q = CV% an d B = I^ 1 — 1 — a. Equation (|A8|I has 
real solutions for both B < 1 (then Sf < 0) and B > 1 
(then 5/ > 0). Reduction of the pulse width in the first 
case corresponds to the converging laser pulse, while in 
the second case the laser diverges. The border of laser 
pulse stability against the transverse distortions, B = 1, 
gives the implicit equation for the stability point, 



I(£ c ) = (2 + a)- 



(A9) 



The point £ c separates diverging and converging parts 
in the pulse profile. The function /(£) grows monotoni- 
cally from the leading front towards the trailing edge of 
the pulse. Therefore, the pulse portion which lays be- 
tween the leading front and the point £ c spreads, while 
the part beyond £ c focuses. As a result, in the pro- 
cess of propagation the short pulse acquires the form 
similar to "beet-root". In the vicinity of the thresh- 
old point £ c , the pulse width given by Eq. (|A8|I reads 

We examine the case of a Gaussian laser pulse, Aq (£) = 
a exp[— 2 ln2£ 2 / (ct l ) 2 }, with the center located at £ — 0. 
Substituting Aq(£) into Eq. <|A6|> gives 



notations are the same as in Eq. QJ. At the initial mo- 
ment t = 0, the pulse center £- = resides at z — —zq. 
The pulse propagates towards positive z. 

In the linear approximation, the frequency of electron 
plasma oscillation equals uj p , and the phase velocity v p of 
the wake coincides with the group velocity v g of the pulse. 
In the weakly nonlinear approximation, the frequency of 
plasma oscillations is downshifted because of relativistic 
increase of the mass of oscillating electron 0] . The in- 
tensity of the Gaussian beam I|B1I) varies along the laser 
path, and so does the amplitude of the plasma wake. 
Hence, relativistic reduction of the plasma frequency al- 
ters as a function of propagation distance, and this phase 
slippage characterizes the variation of the phase velocity 
of the wake across the plasma. Consideration of this ef- 
fect is easier in the case of a wide, k p ro 3> 1, and not so 
intense, clq < 1, laser pulse when the wake electric field 
is the potential one. Then, for the laser pulse amplitude 
taken in the form I|B1|) , the dimcnsionless wake potential 
4> = etp/(mc 2 ) is given by Q 



(a /2) 2 g(x) 



1 + C 2 



(B2) 



x exp 



2r 2 



r 2 (l + C 2 ) 



Aoj p (z, r) 



where 



8 In 2 



where g{x) = Xy/ir/2 exp(— x 2 /8) depends on the dimen- 
sionless pulse length x — ui p tl /V2\n2, the relativistic 
frequency shift pjj is Auj p = -(3/16)cj p [(a /2) 2 /(l + 
dx[l-&(x)], (A10) C 2 )] 2 exph(2r/r )7(l + C 2 )].9 2 (2;). The on-axis phase of 



-&Vln2/(cT L ) 



the wake potential follows from Eq. (|B2 



2 f x 

$(.t) = —= \ dtexp(-t 2 
V n Jo 



kp(v g t - z) I 1 



3 

Tg 



(a /2) 2 g(x, 



n 2 " 



l + C 2 



is the probability integral. Substantial modification of 
the laser pulse shape occurs when the threshold condi- 
tion IjAlOJ) is met at the pulse center (i.e., £ c = 0). In 
consistence with Eq. ijAlOfl . this criterion takes the form 



Then, the wake phase velocity v p reads 



36/ dt 

"befdz 



(ao/2) 2 g(z) 

l + C 2 



a/z R 

l + C 2 



(B3) 



(B4) 



41n2 



0()c 



u p t l V 1 + (fc p r /4) 2 



(All) 



APPENDIX B: THE WAKEFIELD OF AN 
ULTRA-SHORT FOCUSING LASER PULSE 

The wakefield is excited by the ponderomotive force 
associated with a laser pulse. Under the assumption that 
the linear theory of diffraction of Gaussian beams [l2T ] 
holds, the laser pulse envelope a(r, reads 



Of) 



: exp 



-21n2- 



e 



(cr L y 



r 2 {V 



-c 2 ) 



(Bl) 



The quantity t;/v g characterizes a positive time delay 
with respect to the pulse center, and C characterizes the 
laser pulse position relative to the focal plane. When 
the pulse passes the focal plane, C becomes positive. 
The wake phase velocity coincides with the pulse group 
velocity only far away from the focal plane and ex- 
actly at the focal plane. The wake phase velocity de- 
creases while the pulse approaches the focal plane, and 
reaches the minimum value for z — where 
v p w v g {l - 0.192[(a /2) 2 ff (2:)] 2 C/^}. The wake phase 
velocity exceeds v g and can exceed the vacuum speed of 
light while the pulse moves away from the focal plane and 
radially spreads. Small variation of the phase velocity v p 
with respect to v g can, however, considerably modify the 
Lorentz factor 7 p o = (1 — Vp/c 2 )^ 1 ^ 2 if compared with 
7 S = (1 - v 2 /c 2 )-^ 2 . Eq. CEp gives 



where C = z/Zr, 
C 2 ) — arctanC, v, 



C - v g t -{z + zo), * = (r/r ) 2 C/(l + 
is the group velocity of a pulse, other 



1 P « lg +~ 



(a /2) 2 g(x) 



l + C 2 



i + C : 



2 ' 



(B5) 
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The electrons are injected at negative £ where v p < v g . 
Hence, the resonance with the nonlinear wake field needs 
less electron energy (mc 2r y p o) than the resonance with 
the laser pulse (or a linear wakefield), mc 2 j g . 

Although Eq. (|B5(I is weakly relativistic, we use it 
for a qualitative estimate of resonant gamma-factors for 
different wake periods for the simulation parameters of 
Sec. US and |V| (uj p t l = 0.56, a = 1.72, k p r = 6.3, 



Zb = 4 cm, 7 g = 125). The accelerating phase of the 
second period (See Fig.EJ corresponds to k p £ « 6.5 (i.e., 
£ = —0.01 cm), and, at the point of injection £ = — 1, 
has a gamma- factor 7 p o ~ 80. For the accelerating phase 
of the fifth plasma period (fc p £ » 25) Eq. I|B5|I gives 
7p0 ~ 49. Graphic estimates of resonant Lorentz fac- 
tors in Fig. b) exceed the weakly relativistic estimates 
roughly by a factor two. 
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